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LONGITUDINAL OSCILLATION OF AN AIRPLANE.* 
PART I - PROBLEM AND METHOD. 
By R. FuchB and L. Hopf. 

Introduction . 

All a exo dynamical calculations, wMch form the "basis of the 
designing and analysis of airplanes, are fourjded upon an assu-orp- 
tion of uniform steady rcotion. There are, ho^?ever, many prob- 
lems of practical importance which cannot be disposed of by 
means of such calculations. To these belong all questions re- 
lating to maneuverability and to the maximum stresses undergone 
during flight. Such -oroblems have hitherto been solved by apply-, 
ing established scientific principles relating to steady flight, 
as conformably to the facts as possible and, vt'nexe this method 
Tfas not practicable, by relying in individual oases, on the 
opinion of the pilot. If, however, aerodynamics is to afford a 
v;ider basis for the art of flying, it must elucidate the prob- 
lems of accelerated and disturbed motion, phenomena in an acci- 
dental or intentional disturbance throu^ the deflection of the 

rudder or elevator, or any other change in the conditions of 

flight. 



* From Technische Berichte, Volume III, No. 7, pp. 317-330- 



Our experimental l<noTirledge is sufficiently extensive to af- 
ford, in irany instances, the necessary "basis for the mathematical 
analysis of these phenomena and, where such is not the case, the 
theoretical investigation of the prohlem can indicate the neces- 
sary experiments on models and actual flight. 

The treatment of the whole problem may, at first, appear 
hopeless to the theorist- The problem is that of a body with 
six degrees of freedom moving in a fluid (air) and on Tnrhioh 
forces are acting, whose relation to the position of the body is 
only known empirically- Bryan's great service is the circum- 
vention of these difficulties by the method of small oscilla- 
tions and the opening of a way for the treatment of these prob- 
lems, even though the method is restricted to simple conditions. 
The method of small oscillations is only applicable to condi- 
tions in the neighborhood of a feaown state of equilibrium and, 
as applied to the present instance, is as follows: 

An airplane is in steady flight along a given line. The 
quantities which determine its position and conditions of flight 
(namely, velocity, angle of attack, slope of the flight path, 
angle of bank, rate of side-slipping and curvature of flight 
path) are all interdependent, when the engine generates a defi- 
nite propeller thrust and the rudders have a definite position. 
The values of the above quantities are determined by the condi- 
tions of equilibrium in steady flight, vrhich have never been 
satisfactorily discussed. If one or more of these variables 
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have not the value required by the conditions of equilibrium, 
then the oscillation cannot be steady. It is simply assumed 
that none of these variables differs greatly from the value cor- 
responding to equilibrium in steady flight. 

All terms of the equations of motion, which are due to dis- 
turbances, are -Qien expanded in powers of the small displace- 
ments and only the first term of each series is retained. 
The equations in this form are linear and easily solved. No 
further fundamental difficulties present themselves and only the 
mathematical K)rk (which is som,etimes very hard), has to be per- 
formed. 

All previous calculations refer to small departures from 
rectilinear flight, for which there are two independent groups, 
each consisting of three equations of motion. Longitudinal and 
lateral oscillations take place, in this case, independently of 
one another. Changes in speed, in the slope of the flight-path 
and in the angle of attack produce no lateral or unsymmetrical 
motion oscillations. Banking, side-slipping and yawing, so long 
as they are small, have no effect on the symmetrical oscillations 
which are determined by the three above-mentioned variables. 
On the contrary, it is impossible to separate still further the 
lateral oscillations and treat rotations around the longitudinal 
axis (rolling), separately from rotation around the vertical 
axis Z, (yawing). 

In practically all contributions which deal wiiii the problem 



tlnas simplified, only the question of stability in rectilinear 
flight is raised. The actual course of a disturbance is not 
followed out, but only the question is raised as to whether an 

airplane ipill finally return to equilibrium from the disturbed ^ 

condition due to the small variables introduced, or whether it 
has a tendency to diverge further from it. In the former case, 
the flight condition is termed stable, in the latter, unstable. 
The importance of gind the effect due to aerodynamical quajatities, 
as determined by the design of the airplane, are brought out by 
this procedure, but, on the contrary, it is not cDear as to v/hat 
significance is given the term stability and what qualities the 
airplane will exhibit in the hands of the pilot. Keither can 
it be maintained that the aerodynamic theory of stability has 
borne fruit in practice, except possibly in England, where it 
has been supplemented by systematic tests on models. So long as 
the stability of only one flight condition is examined, all the 
above-mentioned questions remain undecided. Above ail, we do 
not understand how stability and majtieuverability are mutually 
related and whether an airplane can be handled as well when con- 
structed with a high degree of stability, and how an airplane re- 
acts to an accidental disturbance. 

Reissner* first recognized the need of investigating aore 
deeply into the actual facts and of going beyond the question of _ 
mere possession of stability to describe the actual course of the 

disturbed motion. Reissner and his pupil Gehlen** were the first 

*Zeitschrift fiir Fluetechnik und Motorluftschif fahrt, 1910, 
Nos- 9 and 10. — 
** "Dissertation," r)ublished by R. Oldehbourg, Munich, 1913. 



to work out the problem for lateral oscillations. Gehlen not on- 
ly solves the equation of stability -when longitudinal and lateral 
oscillations are mtually indepencLent, but also determines the 
integration constants, from arbitrary initial, values of a dis- 
turbance or of a Eovenient ox a rudder, and gives a complete de- 
scription of the phenomena involved, namely: the angle of bank, 
the radius of the turn, and the drift of the airplane. 

This analysis, which represents the utmost that can be at- 
tained by the method of small oscillations, still has the defect 
that it is restricted to small deviations from rectilinear flight. 
In curve-flight the problem is not one of small deviations -ohly. 
The problems cannot be solved without taking into account, at ' 
the same time, the longitudinal oscillations, yrhioh are separated 
from the lateral oscillations in the method used for small oscil- 
lations. For instance, the question today depends on vrhether an 
airplane is ascending or descending in a turn. For finite oscil- 
lations, which may differ to any extent from rectilinear flight, 
even now the problem cannot be attacked, since not even the gen- 
eral case of steady motion which includes both rectilinear and 
curving flight, has yet been solved. 

The question of longitudinal oscillation is different, since 
there are, of course, steady longitudinal, oscillations v/ith any 
desired velocity, angle of attack ana. slope of flight path, with- 
out lateral oscillations. Such a steady forward motion presents 
many points of technical importance, some of which we will touch 
upon. There is, for instance, the question of static stability. 



on which, something has already been published in the Technische 
Berichte.* It does not. follow, from static considerations, as 
to v;hat magnitude of static stability should be chosen for an 
airplane, according to its purpose, and e;ven analysis by the 
method of small oscillations furnishes no conclusion on this 
point. It is knoTim that airplanes which are ordinarily unstable, 
may nevertheless fly well under special circumstances. It is 
also known that even with a stable airplane, conditions may arise 
out of which the airplgucie can only be rescued with difficulty* 
The best Isnown example is "stalling," in which the aisrplane no 
longer obeys the elevator. Although the result is usiially an 
unsymmetrioal oscillation, it is generally begun by "stalling" 
(that is, a symmetrical condition of flight into which enter 
hitherto unexplained relations). Fnenomena have also been ob- 
served in diving, which may give rise to danger- The problems 
of stressing are purely problems of longitudinal oscillation, 
since airplanes have hitherto only been tested under symmetrical 
loads. We have to determine what centrifugal forces appear in 
flattening out after a dive, or what lift coefficients and ve- 
locities combine, in passing from high speed at small angles of 
attack to low speed at large angles of attack. 

There are, as yet, no analytical methods for longitudinal 
oscillations as valuable as those of Gehlen for lateral oscilla- 

* Technische Berichte, Volume I, No. 1, p. 16, and following; 
Volume I, Ko.4, p. 108 and following; Volume II, No.l, p. 33; Vol- 
ume II, No. 3, p*463 and following. 
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tions.* Instead of extendliig these calculations, based oti the 
method of small oscillations, the range of which is difficult to 
perceive, and which lead to most complioated calculations, we 
have adopted another method which is not restricted to small os- 
cillations. This was encouraged partly by the results of a num- 
ber of purely practical numerical computations (according to . 
T:hich, approved airplanes appear to have almost neutral equilib- 
rium) and partly by Lanohester's "Theory of Phugoids," which 
includes general longitudinal oscillations differing widely from 
steady fli^t. 

Lift is considered as the only air force in the phugoid theo 
ry. Drag is neglected, thus eliminating all dissipative forces, 
and the principle of energy supplies a simple solution of the » 
equations of motion. The angle of attack is further assumed to 
be invariable during the phole'period of oscillation. In this 
rray, all empirical relations are excluded and the whole motion 
may be analyiiically presented. The significance of these simpli- 
fications Tfill be gone into in a subsequent section. However 
bold they may seem, the result nevertheless agrees with motions 
actually observed in flight. "Looping the loop"' was recognized 
in the phugoid theory long before Pegoud- Paper darts, such as 
children play with, and gliding models, thrown into the air, de- 
scribe Eotions which agree exactly with those required by the 
phugoid theory. Actual airplanes, when left to themselves, do 
not, however, fly in phugoids and the suppositions of the phugoid 

theory mu st, therefore, fail entirely in the domain of full-size 
* Papers by Bryan and his pupils, ^/hich are difficult to obtain, 
appear, however, to deal with this subject. 
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controllable aircraft. 

It has been pointed out by Von Karican and Treff tz* tliat tlie 
phugoid oscillation gives a solution for the ordinary equations 
of longifudinal oscillations when the static stability is infinite. 
An airplane then resists every change in its angle of attack and 
the most iiirportant aseuEsption of the phugoid theory is satisfied- 
The possibility of controlling the airplane disappears completely^ 
since transition from one condition of fligjht to another is incon- 
ceivable, without changing the angle of attac'^i. , Since no airplane 
can be built TThich vrt.ll describe a phugoid, the phugoid theory 
is not suitable for the elucidation of all these relations- It 
fails to solve the most important problem of all, namely, th^t of 
controllability. Ifuraerioal calculations, in fact, lead to the 
anticipated result that the equilibrium of serviceable airplanes 
is not infinitely stable, but, on the contrary, is very small 
(positive or negative) and that airplanes, in the first approxi- 
mation, are nei^-^ral (Technische Berichte, Volume II, No. 3, p. 463) 

As will be shown later, neutral equilibrium greatly simpli- 
fies the equations of oscillation. Ib the simplest case, indeed, 
it is not the angle of attack that remains unaltered, but the 
angle between the longitudinal aris of the airplane and the hori- 
zontal. The other equations can be easily solved, but not in a 
closed expression, as in the phugoid theory, since the forces only 

dep end empirically on the noTT variable angle of attack. The solu 

* dber Langsstabilita t un6. Langs schtringungen von Flugzeugen" 
(Longitudinal stability and longitudinal oscil^-ations of airplanes) 
Jahrbuch der Wi s e en sohaft lichen Gesellschaft fur Luftfahrt, Volume 
III, 19i4-15, p. 116- 



tion, however, exhibits one very definite charaot eristic, namely, 
the oscillation ■which it represents hae? trro distinct phsses with 
respect to time, in that the forces at right angles to the flight- 
path first reach eqailibriun and then (imch nore slowly) , the 
foi'ces in the direction of flight. 

Our further considerations are based on this fact. The re- 
sult is used suggestively, in an attempt to. find an appro xiira te 

solution of the non-neutial airplane having the above-described 
characteristic. It is assumed that the velocity changes more _ 
slowly than the other terms determining the conditions of flight 
The method has in every case proved applicable. 1-p leads to a 
step-by-step integration of the equations of oscillation from the^ 
original condition, but the steps are so long that rarely more 
than two are required and, within the range of each step, the in- 
dividual quantities are obtained in closed form, as solutions of 
linear differential eqiiations. 

3y this method, it is easy to represent the most important 
part (the initial stage) of the course of an oscillation without 
being limited to fixed conditions. A variable angle of the con-_ 
trol surface, slow application of the controls or a back smd fortti 
motion of the controls can, in this way, be as easily expressed 
as any accidental external disturbances, gusts, etc In the first 
part, the method will be worked out and the formulas s^ven, while, 
in the following jarts, definite problems will be treated. 
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Pymbols Used in Oscillation Eq-aatlons . 



X = Angle formed with horizontal by tangent to flight-path 
(Fig. 1); 

9 •= Angle "between upper wing and horizontal (Fig- l) . 

a = Angle of attack, angle "between chord of upper wing and 
tangent to flight-path; 

cp = Angle between propeller axis and flight-path 
hence a - cp = i^; ■ 

ir; = Angle of incidence between upper wing and propeller axis, 
% = e - a (Fig- 1); 

W = Dotal weight of airplane in Isilo grams; 

T = Propeller thrust in kilograjEs; 

S = Area of supporting surface in n^? 

V = Resultant speed of airplane in meters per second, con- 
sidered positive in the direction of flight; 

?v = Specific weight of the air in kilograms per ; 

g = Acceleration due to gravity =9.81 meters per second^ 

q = Dynamic pressure = i x — x V , 

2 g 

Ol and Cj^ coefficients of lift and drag; 

L = q S Ox, (lift); D = q S G]3 (drag); 

M = Moment of forces of airplane about its center of gravity. 



measured in sndh. a way that a moment is positive when it 
turns the nose of the airplane downward. The positive 
direction of the turning moment is, therefore, opposite 
to that of the angles 9, cx, Tt- . 
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Forces in the Direction of Flight .- In the direction of flight 
there act: a component of the propeller thrust, T cos qj , a compo- 
nent of the force of gravity, -W sin X ; and the air resistence, 
-Oq q S. We therefore obtain 

]il ^ = T cos® - W sinX - ^ Gd ^ S - . (l) 

g dt ^ 

Forces at Right Angles to the Direc ts nri of Flight-- There act 
at right angles to the direction of fligfit: a component of the 
propeller thrust, T sin cp ; a component of the force of gravity, 
-W cos % ; and the lift, Q S. To this imiet be added, in curve- 
flight, the centrifugal force | ^ in which r is the radius 

g ^ 

of curvature of the flight-path. To determine r, we have .-_ 
(Fig. 2) the equation 

= dX 

r 

The centrifugal force acts in the direction of gravity, ^rhen 
d X is positive. If the direction of gravity is considered neg- 
ative, the centrifugal force is -written 

I V 
" g d t 

Since there are no components of velocity at rigiit angles to the 
flight-path, we obtain 

0 = ^ V ^ + T sin CP - W cos -X + ^ I Ol S V= (H) 

g do 5 ^ 

Moments .- If k is the radius of gyration of the airplane 
about its center of gravity, the positive direction of the moment 
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"being opposite to that of increasing 9, we have 

w die = _ M 
g ^ dt= 

In finding the value of the moment, it Tras borne in mind 

d. 9 

that M depends upo?i V,a and the angular velpoity ^ . 
The resultant moment M, is principally made up* of the moments 
of the wings, together with that of the horizontal stabilizer 
and elevator. Assuming as us\ial, that the moment is proportional 
to the squa.re of the velocity and is a linear function of a, 

1£ = (mo + m^ a) 

and m, V® = ~, the so-called static stability. mQ is deter- 
mined by the position of the elevator at the time. This, in 
turn, fixes the angle of attack at which the moments are in equi- 
librium. It cannot, of course, be assumed that m^ is the same 
for all conditions of flight. If, however, the equations are 
taken for successive intervals of time, within which the angle 
of attack does not vary too much, it may then be safely assumed 
that is constant for the duration of such an interval. An 

interval can, in any case, last only so long as the position of 
the elevator does not vary, that is, so long as mo has the 
same value. 

When 9 itself varies with the time, the moment is a func- 
tion of this variation, ?x = 9 • To obtain the differential 
* dt 

* Zur Berechnung der Langsmomente von Flugzeugen (Calculation of 
tlie longitudinal rraments of airplanes), Teohnische Berichte, 
Volume II, Sro.3, pp. 463-483. 



coefficient of this relation, the damping, it must Tdo xemoaibered 

that any variation of 6 affects a as well as V. The effect 

on V may be neglected, since it is very small, and even the 
alteration of a need only be introduced in the calculation for 
the moment of the horizontal stabilizer and elevator. An estimate 
indicates that the damping factors, which allow for the variation 
of V, only amount to ■— and the damping effect of the wings 
to only the damping of the horizontal tail plane being 

taken as 1. If, therefore, Mg is the moment of the horizontal 
tailplane and elevator, it is only necessary, in the expression 
for K, to add 

dMg ^ ^ dMa da J 
de da de 

In order to be able to express the ratio let r^ indi- 

cate the distance of the middle of the horizontal tail surfaces 
from the center of gravity. By a rotation dG, around the axis, 
the horizontal tail surfaces are lowered by r^ dG . If the air- 
plane simultaneously advances a distance of dx, the angle of 
attack increases by A a, so that 

rrrde rue 
tan A a = = 



dx V • _ 

The cosine of the small angle which the horizontal tailplane 
makes with the propeller axis, is here taken as unity. 7/ith the 
small size of the angle in question, the tangent may be replaced 
by the arc and we then obtain 

da _ rj£ 
de ~ V 



The eqixation of moments thus tates the form 

0 = (f - ma) - d V H (III) 



in vrhich 



= _i_ A M_ d = ^ CD 



Wk da Y Wis: da V® 

while f is a constant determined by the position of the ele- 
vator at the time. 

In equations (l) and (II) we put X = Q -.a; further be- 

g _ 

cause of the small angles: cosa= costp = 1; sin a -■ 5773"^ ana 
sine? = 7=773; all angles being meaisured in degrees. On putting 

^ = 7, the differential equations become 

g=f - gsine + ^oosea - ^ | Oj) S (la) 

II = 7 - (a - iV) + 57.5 g cose + g sinQa]- Gj^S^- 

ilia.) 

^ = 'Y (Ilia) 

^ = (i - ma) - d 7 V (ITa) 
d \t 

k few ad.ditional remarks may here be made on the analyti- 
cal expression for the propeller thrast- If we assume* that 

the thrust decreases as the square of the velocity (which is 

* An exact basis for this law is unforturiately lacking, as yet, _ 
The attempts by Sann (Tcchnische Beridhte, Volume I, IT0.6, pp. 
232-241) VTOuld here be too elaborate. The present assumption 
agrees approximately ^Tith the expression derived by Sverling 
(Zeitschrift f^ir Plugtechnik und Motorluftschiffahrt, 1916, p,127, 
equation ( 8 ) . 
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only approxiaate) or, in other words, that tlie fhrust is repre- 
sented ("witliin the range discussed) "bj the i)arahola — 

T = To - pV^ (3) 

(in ■which, of course, Tq is not exactly the thrust when stand- 
ing, still) then p is still a function of the density of the air. 
We VTill, however, assume that the altitude of the airplane does 
not change materially T7itbin the limits of one of the time in- 
tervals, considered, and that the density of the air therefore 
remains approximately constant* In equation (Ila), T may he 
considered constant, because ^ (a - ivr)V is significant in 
comparison ^ith - ^"^^^ ^ S V and also because the propeller 
thrust plays no part in this equa.tion. In order to adapt equa- 
tion (2) to equation (la), vre tiien write 

T = Oi, - i I S V=P' 

and consider p' combined with Op, that is, v/hen flying with 
the engine going, the coefficient of drag is correspondingly in-, 
creased* T can then be considered invariable in the differen- 
tial equations. Moreover, an increase in Gjj mu.st also be made 
for steep gliding flight, v/hen the propeller is running light. 

II. Airplane Trith Neutral Equilibrium . 

V/hen static stability is zero and the moment does not change 
with the angle of attack, 

da 
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In order that the airplane csay he in equilibrium and ijhat 
the moments set up by the -wings and the tail may balance, it -is 
necessary for the position of the elevators to be chosen so that 
f = 0. The airplane is then in neutral equilibrium and equation 
(iVa) becomes 

t 

d-Y - / d V dt 

|i. = - d7 V, also 7 = G e o (1) 



Assuming that when t = 0, T = ^ = 0, it follows that 



d£ 
dt 

0 = 0, that is, 7=0. In this case, therefore, the angle 
which, the upper wing makes with the horizontal remains unaltered. 

Equations (Ilia) and (iVa) are now eliminated and the vari- 
ations of V and a are determined by 

dt " ^ - S Bin 6^+ gf-g cos %a - 1^ S Od V^, (2) 
da _ 57>5 g cos 6., 57.5 ^- ^ p ('^) 



The terms ^ (a - i^) + g sin a. a of equation (lla) are 
omitted in equation (3), since they are very small in comparison 
with 57.3 g cos 6^ . This can always be done when Qq does not 

approach ±90° . 

From equations (2) and (3), Y and a must now be calculat- 
ed as functions of t. It is convenient, in many cases, to plot 
simultaneous values of a and V as coordinates. The following 
rule is important for the discussion of the resulting V,a curve. 
For each G^, there is, in the V, a plane, a curve, along which 
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there is equilibxium of the forces in the direclion of flight. 

This curve is obtained by putting fr = 0 in equation (2) and 

ax — 

then expressing V as a function of a. For each e^, there is 

also a curve, along which there is equilibrium of the forces at 

right angles to the direction of flight. This curve is derived 
from equation (3) by putting |f = 0- There is coisplete equilib- 
rium of the forces at the point of intersection of these two 
curves. Figs. 3 to 5 show these curves for three different values 
of e^, so that there is equilibrium when a = 3° in Fig. 3, 
a = 10^ in Fig. 4, and a = 15° in Fig. 5. We have taken 
W = 1,530 kg; S = 41.3 m^; t = 485 - 0-05 x ^ x - x s V®; 

' 2 ^ "g" = 1572 ' "^^^ values for Cl and Oj) are taken 
from the polar diagram of the Dfw 0 V. The two equilibrium 
curves intersect in the V,a plane in a second point, in addition 
to the point for which they axe calculated. In the later inves- 
tigations of the so-called "stalled" condition, the importance 
is shown of the question as to whether this second intersection 
is at a greater or a smaller angle of attack than stalling. We 
must, therefore, compute immediately at what value of 6o the two 
curves of equilibrium 3ust touch each other. If we write the 
equations, obtained by putting = 0 and ^ = 0 in equations 
(3) and' (3), in the following form 

we obtain, for the desired point of contact. 
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da da ' ^ • 

From the second of these equations we obtain, for each a, the 
corresponding e^, of complete equilibrium. We further find 

xdOo dBo/cLa da da 
and, consequently, =0- If, therefore, is plotted 

against a for all positions of equilibrium, thi s curve mast have 
a maximum value for the point of contact. 
From equation (4) it folloirs that 

(6) 

da \0i^/ 

The numerical factor 57.3 mst not be forgotten, when a is ex- 
pressed in degrees. 

III. Analytical Calculation of Neutral Equilibrium in Flight . 

In Figs.. 3 to 5, using the data for the airplane Dfw 0 V, 
in addition to the curves of equilibrium for several solutions 
of differential equations, simultaneous values of 7 and a 
are plotted by means of a numerical integration. From all these 



curves, it will be seen that (for a point V, a at some distance 
from the curve ^ ^ °^ equilibrium of forces at right 
angles to the direction of flight) the corresponding integral 
curve always runs almost parallel to the a axis, so that varia- 
tions in a correspond to much smaller variations of Y. In 
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other words, the forces at right angles to the direction of flig ^^t 
attain a state of equilibrium much more rapidly than those in the 
direction of fli^t. The reason is that much smaller values ap- 
pear on the right side of equation 2 (II) than of equation 3 (III). 

Only when we approach the curve of equilibrium Ji= = 0. does the 

at ' 

order of magnitude of ^ and || become the same. 

Accordingly, in the analytical calculation of the variations 
in the velocity and in the angle of attack, each step will be di- 
vided into separate intervals of time. 

Case A.- Let us suppose that the equilibrium of an airplane 
has been disturbed, so that, at the beginning of the unsteady 
flight thus initiated, the velocity Vq and the angle of attack 



% determine a point far removed from the V,a curve = 0. 

The return toward neutral equilibrium, must first of all be exam- 

ined, as to the time when the forces at right angles to the line 

of flight are appro :i:imately in equilibrium. As a first approxi- 
mation, we put V = and then determine, from equation 3 (ll) 

da ^ 57.3 g cosGj, _ 57^ q tt r (l.) 

In carrying out this integration, only the time intervals 
will be considered, during which G-^ may be regarded as a lin- 
ear function of a . This is possible with G-^ through a wide 
range, unless we come quite near the maximum lift. In the neigh- 
borhood of this maximum value the rangg becomes smaller and the 
time intervals of the integration must be taken correspondingly 
smaller. 
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If we pat for Cj, in equation (l) ■ 

Cl = CLo + ClLia (2) 
and further, for shortness 

c = 5p I (3) 
equation (1) then takes the form 



= g - e.CLo Vo - V^a (la) 



da 



dt Vq 



The solution of this eqtiation, which has the value 
for t"= 0, is 

a = L + (a, - De"^^^^ (4) 



o 

when 



^ _ 57.3 cos _ Clq. 

In order to make a second approximation for V, this value 
is put for a, together with V = Tq* on "^^h® right-hand side 
of equation 2 (II) . 

^ = ^ - g sin e + a - Gd (5) 

dt W ^ 0 57.3 57.3 ^ o 

and takes, for simplicity in the time interval under considera-. 
tion, as a linear function of 

A less simple expression offers no difficulty, but makes 
,the result less concise,, without affecting it materially. This _ 
gives 
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dt ^? S ^^^^o^ 57-i3 : 57.3 57.3 



We find, by integration, the following sol-ution which has 
the value V = Vq for t = 0: 

V = V - p + Nt + P e"'°^^ (8) 

when 

p = ( ^0 _ -g ooseo > (a^ - L) 

^ 57.3 ClI 57.3 €0^1 Vq / ° 

N = ^ - g sine '^QBgo _ e GDI V^^ N ^ c_^a_Ial . 

^ ^'o V 57,3 57.3 57.3 

* 

Example : In the ahove numerical example, S = 41.3 m^; 

W = 1530 kg; T = 485 - 0.05 x ^ x - S V^, 

S 

To the values of Cj) in the polar diagram (Dfw 0 V) , there mast 
be added, in accordance with the rule for T, the amount 0,05 
and also the coefficient of structural drag 0.0336. We may then, 
according to the dimensions of the model, put 

Olo = 0.325, = 0,0673; = 0.115, Cdi =.0.00562 

For 0^= 7°^ there is equilibrium when a = 3°. 
T/e obtain: 

a = 103m _ 4.84 ^ (a^ _ 103gO ^ ^.g^N ^-0.05*3 Vpt 
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V = Vo - (0.00146 Vo - f a„ - + 4,84^ + 

+ (0.27 + - 0*00124 V'^Jt + (10) 



+ (0.00146 V- - (.a - 10300 ^.84^ e 



-0.06 43 VQt 



Vo ^ " Vo 

These expxesaionB will be discussed later in a mmerioal 
example . 

Case B.- If the point determined by Vq and Cq in the _ 
V, a. plane is very close to the V, a {^'- 0) curve at the coffl- 
mencement of the motion, we then find, in contrast to Case A, that 
the variation of Y is of the same order as that of a. and we 
must consider the two equations together 

dV To* « fi" 

dt = If ~ S sin + cos e„ a - 

(ODo + Cdi a ) = X (V, a), ' (11) 

If = 57.3 s oose ^ _ , (0^^ + c^^cc) V = s(V,a) (12^^ 
If a solution commences close to the equilibrium curve 

do. 

^ =^ 0, as assumed here, simultaneous values of V, a will remain 
close to it throughout. To make this clear, we will consider the 
solution again in time intervals, within which X (V,a) and 
b(V, a) may be regarded ss linear (which is obviously always pes- 



sible) and we will suppose the final values V, a of any time in- 
terval, to be the initial values Vo^a^ of the succeeding inter- 
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val. Examples sHow that, in practice, only one or two time intex- 
vals are required. 

For a single time interval, we have 



da 
dt 



At the beginning of the first interval, s{^o,%) is 
almost 0 and we have 

p. = = - ^ (°Do + %i c^) Vo, 

O 

/ dt\ _ g cos Gq _ ^Opi Vq^ . 
" V dS^/- " 57.3 57.3 

/ds\ 57.3 g cos Qq 



(11a) 
(12a) 



c 

a 



, = =. - € ClI TTq. 

o 



Equations (11a) and (12a) ere solved by putting 



and we find ^ 

L =X (Vo,ao) p^q^!4q^ .B= ^^o, «o q^-p^q- 

r, and r, are the roots of the equation 

- r (p^ + q^) + P, " P3 = ^ 
The factors 0^, O2 , d, , d^ follow from 



(15) 



(16) 
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+ d-s = - B, rgtJg = Qgda + PgCL, , Xgda = iqiCs+p^dg . (18) 

Equation (16) has real negative roots, so long as the initial 
velocity does not fall to stalling speed. In the former example, 
when Yq, lies close to X (Vq,^^) = 0, we have 

r- 4- 5.68 y lO-X + l^Q^ ooBB^ ^ ^ ^ 3^ ^ ^ 

^o 

1.93 X 10^ cos^ 

. V 

and the rcots become complex, only when falls below 

23.4ycos . 

The oscillations from the instant we approach the line of 
equilibrium of forces normal to the direction of flight, can now 
be surveyed in detail. If we again consider siimiltaneous values 
of V, a then the V, a curve can only reach .the line of equilib- 
rium for the perpendicular forces, when V and a rise or fall 

together, since, at the instant of crossing the line of equilib- 

da " 
^i^ni = 0, the curve runs in the direction of the V-axis. 

Should a, for instance, rise and Y fall, then ^ would first 

be positive, then zero and then again positive, that is ^ must 

have a minimm value and, at the same time, become zero at the 

instant when the line of equilibrium is reached- The expressions 



II = c.r.e''^^ H- c,tJ'^ and ^ = cJ, r,^ e''^^ r/ e^^ ^ cannot 

dt 

be zero together, when and have different values- If, 
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therefore, a increases and V decreaseB, or a decreases and 

V increases, then the Y, a curve will certainly remain perma- 
nently close to the line of equilihriun; but -will only reach it 
after a very long time (theoretically t = «) - When a and 

Y increase or decrease simultaneously, the line of equilihrium 
T7ill be crossed once and thereafter the V,a curve will again 
rema,in in close proximity to the line of equilibrium- 

The iTiOtion can be understood better from a numerical example. 
If W 1530 kg, S = 41,3 

0^ = 0.335 + 0.0673a , 
Ojj = 0.115 + 0.00562 a 

then equilibrium exists when a= S°, 9 =7° and V = 36.2 m 

per second. 

Let the equilibrium be so disturbed that, at the. beginning 
of the unsteady flight, % = 5°, = 43.3 m per second (Fig. 6}_ 
are values which give points lying far frora the line of equilib- 
rium of the normal forces. In the first part of the oscillation, 
the solution therefore corresponds to Case A: 

~ 3' 3 5 t 

a = 0.67 + 4.33e (19) 

V = 43.2 - l.llt + 0.0407e"'"^^ (30) 

The calculated values of V and a have been plotted in 
Fig. 6 as functions of each other and in Fig. 7, singly, as 
functions of t. The result is, moreover, compared wi-fii a very 
careful numerical integration (dash lines) and the excellent 
agreement between the curves shows that the analytical calcula- 
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tion is very exact. 

The equilibrium curve is reached in 1.12 seconds, when 

V = 42 m per second and a = 1*^. Now take case B and the vs-lues 
become 

a = 2.65 - 1.77e-°-""^*-'-''^ (21) 

+ o-iae-'-''^^"-'' 

V= 36.05 + 5.95e-°''""^"^^-'""'^ (22) 

Equations (21) and (22) can be used for the whole course up 
to t = <», since, for t = », they give ct = 2.65° and 

V = 36.05 m per second, which, therefore, come very close to the 
coordinate values » = 3° and V =■ 36.2 m per second. In order 
to estimate the time it actually takes to' restore equilibrium, it 
Kust be borne in mind that the term containing ^-^•^''(t-i- is)^ 
(which from the first is vanishingly small in equation (22) and _ 
therefore can be entirely omitted) diminishes rapidly. The term 

_ -0.1S2(t-l' 12) 

has the value 0.1 after 28 seconds. It may, 
therefore, be said that, nith the given disturbance, equilibrium 
is practicslly reached in about half a minate 

IV. The General Case. Discuseion of the Constants. Analytical 

Treatment . 

The general equations for the velocity V, the angle of 
attack -a, the angle of inclination of the upper wing to the hor- 
izontal 6, and the angular velocity " ^ are 
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e = 5p ^ s, (la) 



da 
dt 



- 1 ("is (a - i„) + 57.3 g cos 6 



+ g sin ea] - t Cl V, (Ha) 



= •> (Ilia) 

dt 

(f - ma)V^ - d 7 V. (IV^a) 

Among the coefficients appearing in these equations, some 

are always invariable even under different flight conditions 

(permanent constants), while others vary under different flight 

conditions and can only "be regarded as constant within a given 
time inteinrsLl (temporary constants). 

To the former class "belbng: 

1. Total weight of airplane, neglecting variation in weight 
due, for instance-, to consumption of fuel; 

2. Supporting surface; 

3. Angle between upper wing and propeller axis; 

4. Acceleration due to gravity, g; 

g _ d Mti 

5- Damping coefficient, d = -^H — ^2 
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The temporary constants are: 

1. Propeller .tlirust and air density. 

For the propeller thrust, the expression 

= To - -J- - S p' 
was introduced. T, X and p' are constant within a given time in- 
terval, "but all three quantities may vary in different time inter 
vals. 

2. The static stability m = — is regarded as con- 

Wk V da 

stant within a given time interval, hut it is possible to use dif 
ferent values for m in different time intervals, ^vhen passing 
from one state of flight to another. 

3. As already stated, f varies with the position of the 
elevator. If, for instance, at the beginning of the oscillation, 
the moments of the wings ai:^ of the tail balance at an angle of 
attack of 3*^ and the elevator is then turned so they balance at 
9°, we will have f = 3 m for the first interval and f = 9 m 
for the second. 

4. The coefficients 0^^ and Cp are here introduced as lin- 
ear functions of a-: Cj. = 0-^ + 0-^^ a, Oq = Cjjq + Cjjj^ a, in 
this connection G-^o> ^Ll» ^Bo' ^1 assiuned to be constant 
within any given time interval. These coefficients will, of — 
course, vary in the different time intervals, if a. increases or 
decreases. 
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lf WG now undertake the solution of the general equations hy 
namerical processes, with given values of the permanent and tempo- 
rary constants, starting from a definite instant, t. = 0, with 

it 

arbitrary initial values of Vq, a^, 9^^ '•/Trill always be found 
that at first the velocity only changes slowly, in comparison 
with the angles. This fact offers a very easy way for analytical 
treatment, by assTiraing in the first approximation, as in the case 
of neutral equilibrium, that V = Vq constant- Equation (la) 
drops out and we have only to solve equations (Ila), (Illa), and 
(iVa), in which Y is put equal to Vq. These equations are, 
however, all linear, when 0-^ and Gj-, are linear functions of a, 
which is a great advantage in TOrMng out the problem. Moreover, 
these three equations are reducible to two, provided a certain 
correction is introduced for diving. The values found for 
"^O' °'» ®» P"* ^^""^o equation (la) and we obtain, by 

simple integration, a second approximation for V, which, to- 
gether with the previous values for a, 8 and 7, presents an 
excellent solution for a definite time interval- These analytical 
expressions, as shown by comparison with solutions by means of 
fixed coefficients, give the actual path of flight very well for 
an interval of about two seconds. If it is desired -to follow, 
during an extended period, until equilibrium is reached, the non- 
steady flight caused by any disturbance on a stable airplane, by 
the same methods as for a neutral airplane, the above calculation 
can be used in conjunction with the method of small oscillations. 
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During the first seconds, we calculate by the above method and 
thereby determine -the course of V, a, 6^ 7. When a has com- 
pleted its large variations and has substantially reached its 
equilibriTim value, we calculate the further course, up to equi- 
librium, by the method of small oscillations, making use of all 
four equations and proceeding from the final values of the present 
method. Such an example is worked out in No. VII. 

V- Problems of the General Oase . — 

1. Let an airplane be in equilibrium, with all permanent and 
temporary constants known, and let the equilibrium be disturbed 
by some cause, such as a gust, so that the velocity is changed to. 
Vq, the angle of. attack to Oq, the inclination with the hori- 
zontal to 9^^ and the angular velocity to 7^. What is the 
course of the non-steady motion now set up? More especially, how 
does a stable airplane return to equilibrium? 

3. Let an airplane b.e in equilibrium and a deflection be im- 
parted to the elevator. The values in the state of equilibrium 

a, 6,7 =0 are to be taken as initial values. In the differ- 
ential equations, hotrever, a value of f is to be put, correspond 
ing T7ith the new position of the elevator- Again a non-steady mo- 
tion sets in, whidi has to be followed. It must be especially in- 
vestigated, as to how this deflection of the elevator affects 
stable (m > 0), neutral (m = O) and unstable (m < O) airplane 
By the method of subdivision into time intervals, it is always 
possible to give successively different elevator settings. At the 
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moment of setting the elevator, we must start from the initial 
values of V, a, 9, Y, but, on the other hand, we must introduce, 
into the differential equations, the particular value of f which 
corresponds to the new elevator setting. 

3. Let an airplane be in flight under engine power- At a given _ 
moment, the engine is shut off and a new setting is simultaneously 
given the elevator. The values for V, a, 9 and 7, during engine- _ 
driven flight, stand as initial values. The differential equa- 
tions must, however, be those of gliding flight; that is, T = 0, 

and in the expression = Oqq + Oqio- , for the corresponding 

a position, p', defined by T = J- ^ S V^p' must be omitted and, 

in its place, an amount put, which corresponds to the drag of the 

propeller revolving slowly. The reverse takes place when passing 

from gliding flight to power f light - 

4. The stalled condition C8.n very well be treated by the 
present method, since, precisely in this condition, the velocity 
changes very slowly. All the phenomena peculiar to stalling can, 

therefore, be represented by the general formulas given below, by 

putting the initial values characteristic of this condition (large 
angle of attack and low speed) in the differential equations for 
such values of Cl = Olq + ^LX°-> = <^Do + correspond to 

the a position. 

5- This method also suffices admirably for the treatment of 
diving fli^t, since the velocity in this case has been found to 
change but slow^ly with variations in the coefficients, variation 
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of f due to change of tli§ elevator position and variation of 
T in passing from engine-driven flight to gliding flight. _ 

In all these cases, we have to deal with the following 
mathematical problem. A system of four differential equations 
is given, with definite values, of the permanent and temporary 
constants. Solutions for the four variables are sought as func- 
tions of the time with initial values of Yq, clq, 6o» "Vo*- ^^■^ 
t = 0. If these solutions are to hold for a fairly long period, 
the same problems must be solved for consecutive time intervals, 
the final values of V, a, 6, Y, of the one time interval being 
the initial values of the next. The general analytical method 
of working put these problems is given below and explained by 
examples. The actual solution of the above special problems is 
reserved. 

VI. AioTolication of the Analytical Process to the 

General Problem . 

In all non-steady flights, the velocity V varies but slow- 
ly in comparison with the angles ct and 6. In solving eo.uations 
(la), (Ila), (Illa), and (iVa), it is, therefore, assumed, in the 
first approximation, that V = Vq is constant. We then have 
to deal with the following equations: 

II = 57.3^g cos e _ + ^Ll °^)Vo + 'V ■ (23) 

(24) 

5i = (f - ma)Vo^ - d Yo Y (25) 

at 

in which r, % 

2= As (26) 
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Equation (la) comes first into consideration in seeking a second 
approximation for V. At first, the terms 

- (a - i^) + g sin 8a are again neglected in comparison with 
57.3 g cos e ■• If, for instance, T = 485 kg, TiT = 1530 kg, and 
e = 30° then 

- (a - i^) + g sin ea = 0.34a + 3. 11 i^; 
57.3 g cos e a = 530 a ; 

and this neglect is, therefore, justified. When B approaches 
-90° in diving flight, then these neglected terms again come into 
consideration. The correction, which then becomes necessary, 
will be discussed later. 

It further appears that, except in diving with large oscilla- 
tions of e, 57.3 cos 8 can be replaced by ST'.S cos 8q, in 
which Qq is the initial value of 8, since (with the variations 
in 9 considered here) 57.3 cos 9 only changes by a small per- 
centage, which (as comparison with numerous exact calculations 
always reaffirms) does not materially affect the result. We have, 
therefore, only two equations to deal with: 

If = 57.3 g^cos Q o - c V, - € Cu Vo a + 7 (37) 

^ = (f - ma) - d Vo'Y (28) 

from which 9 is at once given by 

11=7 (39) 
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The solution of equations (37) and (28) is obtained by the 

values ^ . _ . 

a = L + p.e''^^ + Pfe©^ 

7 = B + a e^^ * + q^e^^"^ 

On patting 

W = m + € d Oli (30) 

we find 

L - [ 57.3 g cos eo - ^ ^Lo ^o'' ^ + I ^^^^ 

B = - ^ [57. 3 g cos 00 - e V 3 + ^ ^/^"^ ^^2) 
The values of P, , ISg , , qg , , , are determined by 

Oq - L = iflj + Big ; p, = - € 0^1 pj + q^ ; 
Pgi = - e Oli Vq + q^; 

■Yq - B = nj + n^; q,rj = - m Vq® Pj - d q, ; 
qsTs = - m Vq^ P2 - d Vo~ Qs- 

We therefore obtain for r, and Xg the quadratic equation 

r^ + r Vo(d + e Oli) + N V©^ = 0, (33) 

from which we obtain 



= - ^2 ^o + ^oT^J^* 



Is = - "^o - Vq/r,^ - m (34) 

when T T , V / x 

= I - d), Rs = I (€ Cli + d) (35) 

Hence 
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cx = L -f- e"^=^°^ [a, cos V^, t/ m - i- 

Vot/ m- ], (36) 



+ * sm 



7 = B + e"^^^°^ [e, cos Y^t J m - Rj= + 
+ sin V^t / m - R,0 , 

e = 8o - Cj + Bt + e"^=^'°* cos V^ty^m - Rj^ + 

Vot y m - R/ ] , 



+ i 



sin 



y ra - Ra 
%=e-a=eo-<lj - L + Bt + 

sin Vo't y m - R,^ j 
We have, at the same time, 
= Oo - ag = - Ri (ao - L) + ^ .C^o - B) , 



(57) 



(38) 



(39) 



(40) 



Si = - B, ea = - m Vq (a^ - L) + R, CYq - B), (41) 

= ^[ m Vo(do - L) - € c^i O'o - > . (^2) 

©2 = m Vo RsCtto _ L) + (m - £ Oli R, ) C^o - B)] , 

= - Vo d (a^ - L) + - B] , (43) 

" N Vo 

= -L^rVoCm + d R, ) (ao - L) - R^ CVq - B)] • 

. N Vn 



- 36 - 

These formulas remain unchanged, when m - R,^ < 0, excep 
ing that it is necessary to replace m - R^^ by ^R,^ m and 
the trigonometrical functions "b^"" the corresponding hyperholic 
functions. 

In the case N = 0, ishich is not specially notable in its 
charrxteri sties, the forrmilas break down. If, in such an in- 
stance, we take 

K = [57.3 g cos eo - ^ Clo V 3 + ^^^^ 

= - ^ GlI ^ + € Gx,i K - f ^ 
2 Ra Vq 3 R;^ 

they then become 

a=ao+KVt+A(l -e'^^='^°*), (45) 
o 

— 2 RrjV t - 

y = + e Cj..j K V^^t - Ad Yq (1 - e ' ° (46) 
e = Go + (*Yo - Ad Vo)t +1 £0^1 K t^ + 

+ Ad^'^d - e-~^^o''^), (47) 
2R3 

X .= eo - °-o + CYo - Ad Vo - E Vo) t + 

+ i- £ Gli K Vo^I? - -^^^ (1 - e-^^^^^^"^)- (48) 

A second approximation for V is obtained from 

g Go cos 60 



V = Vq + / dt - g sin Go + 



57-3 



g sin Go ciq Po ^ ^Do "^o^ . / o . o tr 2^ 
57 5^ - 5773 cos Go - c 0^1 ) - 

R c J- S sin Bo do 1 / 
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in which the values found above for a and 0 must "be inserted. 

For 3'udging the course of flight, the angle X , between 
the tangent to the path and the horizontal, is of special impor- 
tance. If we wish to know how the course of fli^t is influ- 
enced by the setting of the elevator, we mst consider ^ , 
that is, the variation of X with respect to f , the variables 
which fix the elevator setting. For this Tare find 



[- 

/ — . Ri^ + - ^ 



2 Rg cos V^t / m - + 



sin V^t ym - R,^ j 



(50) 



or, when N = 0, 



If = C 1 + 2 V„= t= - 3 R, T„t - (51) 

VII. Examples . 

Taking the same data as above, S = 41.3 m^, W = 1530 kg, 
T = 485, -T^ = 0.81 (at an altitude of 3000 m) , 

2 ^ ISTS' ^Lo = 0-325, Cl.i = 0.0672; = 0.115, 

Oj3i = 0.00562; 

the figures correspond approximately to the Dfvr C 7. Let 

m = + 0.00191. (This value was found in calculating the moments 

for the Dfw C V, though with the negative sign) and let 

d = 0.0238 (also the same as for the Dfw 0 V) . 
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1. The state of equilibrium (Vg = 36.2 m per second, 
Og = 3°, 6g = 7° ana. "Yg = .0) is so disturbed that , ths .velocity 
increases to 43-1 m per second and the angle of attack to 6.9°. 
The non-steady motion, which now sets in, is examined and it is 
thus determined in what way the airplane returns to equilibrium. 
In the general formulas, we must put 

= 43.1, = 6.9, = ^g = > '^o = \ = 

Since the moments are assumed to be in equilibrium when 
0^ = 3*^, we must put 

f - 0.00191 X 3 = 0, whence f = 0.00573 
7/e then find, for the first few seconds: 

a = 2.07 + 6.15 e~ ^" ^^cos( 101. 2t + 36-7)°, 

e = 2.69 + 3.24t -5- 4,8e~^' ^^''cos(l01.2t - 27.4)°, 

If = 42.8 - 0.32t - 0.277t^ + 0.3l4e 
cos(l01.2t + 7.8)°. 

After t7.'o seconds, ^e obtain 
a= 1.9°, e = 9°, V = 41.1 m per second, T = 3° per second. 

These va.lues (f rcmp.ining unaltered) can be inserted anew 
in the formula, as . Oq, Sq, V^, 7^, since there is no new set- 
ting of the elevator. We thus obtain for the further course 
( from t = 2) : 

a= 2.27 - 0.5186"'* ^'sin( 224. 7 - 96.5t)°, 
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- 1- e(t- a) 

9 = 9.37 + 3.4(t - 2) - 0.485e 
cos(95.5t - 234.2)°, 

V = 41.1 - 1.03(t - 2) - 0.205(t - sf - 

- 0.0239e cos(96,5t - 172.5). 

The numerical values given by these formulas have been 
plotted in Fig. 8, that is, from t = 0 to t = 2, by the first 
group of forraalas and from t = 2 to t"= 4, by the second group 

The continuous lines have been calculated from the formulas 
while the dotted lines are those obtained by numerical integra-. 
tion. The agreement is excellent. 

In Fig. 9 the same curves (dash) are shovm once- more from 
t = 0 to t = 2; while the course from t = 2 to t = 13 (also 
dash) has been calculated by the method of small oscillations, 
in the neighborhood of the position of equilibrium, with initial 
values corresponding to t = 2. For comparison, the result of 
the numerical integration is also shown (by continuous lines). 

2. Let the airplane be in equilibrium- Then Vg - 36.2 m 
per second, Og = 3°, Qg = 7°, "Yg = 0. Let it be given such 
an elevator setting that the moments are in equilibrium only at 
a = 9°. The non- steady motion set up in this way is to be fol- 
lowed. We again insert m = + 0.00191 and obtain: 

= 36.2, cco = 3^; 6o = 7° 7o = 0, 
f-m9=0, f= 0.0172. 
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We find: 



- 1- *i-t 



cos(85.3t - 43.4)°, 



a = 6.58 - 4.93e 



e = 5.87 + 7t - 3.82e 



- 1K4 1 t 



sin(85.3t - 17.3) , 



V = 36.3 + 0.43t - 0.6t^ - 0.286e 



- 1.-4 1 t 



sin(85.3t - 0.6)° . 



The results of this calculation, from t = 0 to t = 2, axe 
plotted in Fig. 10. 
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Figs. 1,2. 





Figs-1,2. 




Figa.3j^>5' Llnee of eqaU ibrluni of the forces in the direction of flight auad at 
right aniglee thereto . 



Fig. 6,7- 




By analytical intagration 
By numerical integrat-jon 




Time , sec 



Fig. 7 Return to equilibrium of an airplane 
,Tith zero pitching stability. 



Figs. S, 9' 




1 2 

Tiiiie,sec., t 

Fig.S Return to a condition of equilibriun; after 
a distuxbaace during the first four seconds. 



By analytical integration 

By numerical integration 




Fig. 9 Return to equilibrium after a disturbance. 



Fig, 10 




Time, sec. ,t 



igJO Transition from angle of attackja = 3 
to a = 9° produced by a deflection of 
the elevator. 



